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Abstract 

For credit risk management purposes in general, and for allocation of regulatory capital 
by banks in particular (Basel II) , numerical assessments of the credit-worthiness of borrowers 
are indispensable. These assessments are expressed in terms of probabilities of default (PD) 
that should incorporate a certain degree of conservatism in order to reflect the prudential 
risk management style banks are required to apply. In case of credit portfolios that did not 
at all suffer defaults, or very few defaults only over years, the resulting naive zero or close 
to zero estimates would clearly not involve such a sufficient conservatism. As an attempt 
to overcome this issue, we suggest the most prudent estimation principle. This means to 
estimate the PDs by upper confidence bounds while guaranteeing at the same time a PD 
ordering that respects the differences in credit quality indicated by the rating grades. The 
methodology is most easily applied under an assumption of independent default events but 
can be adapted to the case of correlated defaults. 



1 Introduction 

A core input to modern credit risk modeling and managing techniques are probabilities of default 
(PD) per borrower. As such, the accuracy of the PD estimations determines the quality of the 
results of credit risk models. 

One of the obstacles connected with PD estimations can be the low number of defaults, es- 
pecially in the better rating grades. Good rating grades might experience many years without 
any defaults. And even if some defaults occur in a given year, the observed default rates might 
exhibit a high degree of volatility, due to the relatively low number of borrowers in that grade. 
But even entire portfolios with low or no defaults are not uncommon in reality. Examples include 
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portfolios with an overall good quality of borrowers (e.g. sovereign or bank portfolios) as well as 
high- volume-low- number portfolios (e.g. specialized lending). 

Usual bank practices for deriving PD values for such exposures often focus on qualitative 
mapping mechanisms to bank-wide master scales or external ratings. These practices, while 
widespread in the industry, do not entirely satisfy the desire for a statistical foundation of the 
assumed PD values. One may "believe" that the PDs per rating grade appear correct, as well as 
believe that the ordinal ranking and the relative spread between the PDs of two grades is right, 
but information about the absolute PD figures is lacking. Lastly, it could be questioned whether 
these rather q ualitative methods of PD calibration fulfill the minimum requirements set out in 
RCRSI (|2nn4ah. 



The issue has, amongst others, recently been raised in iBBAl (j200J). In that paper, applica- 
tions of causal default models and of exogenous di stribution assumptions on the P Ds across 
the grades have been proposed. In a recent paper, ISchuermann and HansonI (|2Qi]J) present a 
methodology to estimate PD s by means of migration matrices ("duration method", cf. also 
Jafrv and SchuermannL l2004l l . This way, non-zero PDs for high-quality rating grades can be es- 
timated more precisely by counting the borrower migrations through the lower grades to eventual 
default and using Markov chain properties. 

This paper focuses on a different issue of PD estimations in low default portfolios. We present 
a methodology to estimate PDs for portfolios wi thout any defaults, or a v ery low number of 
defaults in the overall portfolio. The proposal by ISchuermann and HansonI does not provide a 
solution for such cases, because the duration method requires a certain number of defaults in at 
least some (usually the low-quality) rating grades. 

For estimating PDs, we use all available quantitative information of the rating system and its 
grades. Moreover, we assume that the ordinal borrower ranking is correct. We do not use any 
additional assumptions or information. Every additional piece of input would be more on the 
assumption side, as the low default property of these portfolios does not provide us with more 
reliable quantitative information. 

Our methodology delivers confidence intervals for the PDs of each rating grade. The PD range 
can be adjusted by the choice of an appropriate confidence level. Moreover, by the most prudent 
estimation principle our methodology yields monotone PD estimates. We look both at the cases 
of uncorrelated and correlated default events, in the latter case under assumptions consistent 
with the Basel risk weight model. 

Moreover, we extend the most prudent estimation by two application variants: First we scale our 
results to overall portfolio central tendencies. Second, we apply our methodology to multi-period 
data and extend our model by time dependencies of the Basel systematic factor. Both variants 
should help to align our principle to realistic data sets and to a range of assumptions that can 
be set according to the specific issues in question when applying our methodology. 
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The paper is structured as follows: The two main concepts underlying the methodology - esti- 
mating PDs as upper confidence bounds and guaranteeing their monotony by the most prudent 
estimation principle - are introduced by two examples that assume independence of the default 
events. The first example deals with a portfolio without any observed defaults. For the second 
example, we modify the first example by assuming that a few defaults have been observed. In 
a further section, we show how the methodology can be modified in order to take into account 
non-zero correlation of default events. This is followed by two sections discussing potential ex- 
tensions of our methodology, in particular the scaling to the overall portfolio central tendency 
and an extension of our model to the multi-period case. The last two sections are devoted to 
discussions of the potential scope of application and of open questions. We conclude with a sum- 
mary of our proposal. In Appendix we provide information on the numerics that is needed 
to implement the estimation approach we suggest. Appendix IbI provides additional numerical 
results to Section El 

2 Example: No Defaults, Assumption of Independence 

The obligors are distributed to rating grades A, B, and C, with frequencies ua, ns, and nc- 
The grade with the highest credit-worthiness is denoted by j4, the grade with the lowest credit- 
worthiness is denoted by C. Neither in A nor in B nor in C any defaults occurred during the 
last observation period. 

We assume that the - still to be estimated - PDs pA of grade A, pB oi grade B^ and pc of grade 
C reflect the decreasing credit- worthiness of the grades, in the sense of the following inequality: 



The inequality implies that we assume the ordinal borrower ranking to be correct. According to 
(|2.1|) . the PD PA of grade A cannot be greater than the PD pc of grade C. As a consequence, the 
most prudent estimate of the value of pA is obtained under the assumption that the probabilities 
PA and PC are equal. Then, from (|2.1j) even follows pA = PB = Pc- Assuming this relation, we 
now proceed in determining a confidence region for pA at confidence level 7. This confidence 
region^ can be described as the set of all admissible values of pA with the property that the 
probability of not observing any default during the observation period is not less than 1 — 7 (for 
instance for 7 = 90%). 

If we have got pA = Pb = Pc-, then the three rating grades A, B, and C do not differ in their 
respective riskiness. Hence we have to deal with a homogeneous sample of size ha + nB + nc 
without any default during the observation period. Assuming unconditional independence of the 
default events, the probability of observing no defaults turns out to be (1 — p^)"-'4+nB+"c_ As a 

^For any value of pA not belonging to this region, the hypothesis that the true PD takes on this value would 
have to be rejected at a type I error level of 1 — 7. 



PA < PB < PC- 




3 



consequence, we have to solve the inequahty 

l-7< (1-^^)"'*+"^+"'' 



(2.2a) 



for PA in order to obtain the confidence region at level 7 for pA as the set of all the values of pA 
such that 

PA<l-(l-7)^/^"-^+"^+"^^. (2.2b) 
If we choose for the sake of illustration 



UA = 100, UB = 400, nc = 300, 



(2.3) 



Table ^ exhibits some values of confidence levels 7 with the corresponding maximum values 
(upper confidence bounds) pA of pA such that 1)2. 2a() is still satisfied. 

Table 1: Upper confidence bound pA ofpA cls a function of the confidence level 
7. No defaults observed, frequencies of obligors in grades given by (|2.3|) . 



7 


50% 


75% 


90% 


95% 


99% 


99.9% 


PA 


0.09% 


0.17% 


0.29% 


0.37% 


0.57% 


0.86% 



According to Table ^ there is a strong dependence of the upper confidence bound pA on the 
confidence level 7. Intuitively, values of 7 smaller than 95% seem more appropriate for estimating 
the PD by pA- 

By inequality (|2.1() . the PD pB of grade B cannot be greater than the PD pc of grade C either. 
Consequently, the most prudent estimate oi pB is obtained by assuming ps = Pc- Assuming 
additional equality with the PD pA of the best grade A would violate the most prudent estimation 
principle, because pA is a lower bound oi pB- If we have got pB = Pc-, then B and C do not 
differ in their respective riskiness and may be considered a homogeneous sample of size UB + nc- 
Therefore, the confidence region at level 7 for is obtained from the inequality 

l-7< (1-Pb)"^+"^. (2.4a) 
(|2.4aj) implies that the confidence region for pB consists of all the values oipB that satisfy 

< 1 - (1 - 7)V("i3+'^c). (2.4b) 

If we take up again the example described by 1)2. 3() . Table [21 exhibits some values of confidence 
levels 7 with the corresponding maximum values (upper confidence bounds) pB oi ps such that 
(tna is still fulfilled. 

For determining the confidence region at level 7 for pc we only make use of the observations in 
grade C because by ()2.1|) there is no obvious upper bound for pc- Hence the confidence region 
at level 7 for pc consists of those values of pc that satisfy the inequality 

l-7< (1-Pc)"^- (2.5a) 
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Table 2: Upper confidence hound pB ofps as a function of the confidence level 
7. No defaults observed, frequencies of obligors in grades given by H2.3|) . 



7 


50% 


75% 


90% 


95% 


99% 


99.9% 


PB 


0.10% 


0.20% 


0.33% 


0.43% 


0.66% 


0.98% 



Equivalently, the confidence region for pc can be described by 

PC < 1 - (1 -7)^/"^- (2.5b) 

Coming back to our example 1)2. 3|) , Table |31 lists some values of confidence levels 7 with the 
corresponding maximum values (upper confidence bounds) pc of pc such that (|2.5al) is still 
fulfilled. 

Table 3: Upper confidence bound pc of pc as a function of the confidence level 
7. No defaults observed, frequencies of obligors in grades given by (|2.3j) . 



7 


50% 


75% 


90% 


95% 


99% 


99.9% 


PC 


0.23% 


0.46% 


0.76% 


0.99% 


1.52% 


2.28% 



Comparison of Tables ^ HI and 01 shows that - besides the confidence level 7 - the applicable 
sample size is a main driver of the upper confidence bound. The smaller the sample size that 
can be made use of, the greater will be the upper confidence bound. This is not an undesirable 
effect because intuitively the credit-worthiness ought to be the better, the greater is the number 
of obligors in a portfolio without any default observation. 

As the results presented so far seem plausible, we suggest to use upper confidence bounds as 
described by (|2.2b|) . (|2.4b|) . and 1)2. 5b|) as estimates for the PDs in portfolios without observed 
defaults. The case of three rating grades we have considered in this section can readily be 
generalized to an arbitrary number of grades. We do not present the details here. 

However, the larger the number of obligors in the entire portfolio is, the more often some defaults 
will occur in some grades at least, even if the general quality of the portfolio is very high. This 
case is not covered by ()2.2b|) . (|2.4b|) . and ()2.5b|) . In the following section, we will show ~ still 
keeping the assumption of independence of the default events - how the most prudent estimation 
methodology can be adapted to the case of a non-zero but still low number of defaults. 

3 Example: Few Defaults, Assumption of Independence 

We consider again the portfolio from Section [2 with the frequencies ha, riB, and nc- In contrast 
to Section [2 this time we assume that during the last period no default was observed in grade 
A, two defaults were observed in grade B, and one default was observed in grade C . 
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As in Section |21 we determine a most prudent confidence region for the PD of A. Again, we 
do so by assuming that the PDs of the three grades are equal. This aUows us to treat the entire 
portfoho as a homogeneous sample of size + nB + riQ- Then the probability of observing not 
more than three defaults is given by the expression 



1=0 



(3.1) 



follows from the fact that the number of defaults in the portfolio is binomially distributed 
as long as the default events are independent. As a consequence of p. If) , the confidence region^ 
at level 7 for is given as the set of all the values of pA that satisfy the inequality 



1-7 < ^('^A+'^s+"c) _^^)nA+nB+nc-i_ 
i=0 



(3.2) 



The tail distribution of a binomial distribution can be expressed in terms of an appropriate beta 
distribution function. Thus, inequality (|3.2)1 can be solved analytically^ for pA- For details, see 
Appendix 1X1 If we assume again that the obligors' numbers per grade are as in ()2.3|) . Table 0] 
shows maximum solutions pA of ()3.2|) for different confidence levels 7. 



Table 4: Upper confidence bound pA of pA o,s a function of the confidence 
level 7. No default observed in grade A, two defaults observed in grade B, one 



default observed 


in grade 


; C, frequencies of obligors in grao 


les given 


7 


50% 


75% 


90% 


95% 


99% 


99.9% 


PA 


0.46% 


0.65% 


0.83% 


0.97% 


1.25% 


1.62% 



Although in grade A no defaults have been observed, the three defaults that occurred during 
the observation period enter the calculation. They effect the upper confidence bounds, which 
are much higher than those in Tabled This is a consequence of the precautionary assumption 
PA = Pb = PC- However, if we alternatively considered grade A alone (by reevaluating (|2.5b|) 
with riA = 100 instead of uq), we would obtain an upper confidence bound 1.38% at level 
7 = 75%. This value is still much higher than the one that has been calculated under the 
precautionary assumption pA = Pb = Pc ~ consequence of the low frequency of obligors in 
grade A in this example. Nevertheless, we see that the methodology described by 1)3. 2() yields 
fairly reasonable results. 

In order to determine the confidence region at level 7 for p^, as in Section [21 we assume that pB 
takes its greatest possible value according to ()2.1() . i.e. that we have pB = Pc- In this situation, 

^We calculate the simple and intuitive exa c t Clap per- Pearson interval. For an overview of this approach, as 
well as potential alternatives, see Brown et al. 1 200lh . 

^Alternatively, solving directly 13.21 for pA by means of numerical tools is not too difficult either (see Ap- 
pendix^ Proposition lA.ll for additional information). 
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we have got a homogeneous portfoho with ns + nc obhgors, PD pB, and three observed defaults. 
In complete analogy to (|3.1|) . the probability of observing no more than three defaults in one 
period then can be written as 



1=0 



riB+nc \ m\ 
i 



(3.3) 



Hence, the confidence region at level 7 for pB turns out to be the set of all the admissible values 
oi Pb which satisfy the inequality 



1-7 < E(""t"^)^'B(i-PBr 

i=0 



+nc—i 



(3.4) 



By analytically or numerically solving 1)3. 4|) for pb ^ with frequencies of obligors in the grades as 
in 1)2. 3() - we obtain Table [S] with some maximum solutions pb oi 1)3. 4|) for different confidence 
levels 7. 



Table 5: Upper confidence hound pB of pB as a function of the confidence 
level 7. No default observed in grade A, two defaults observed in grade B, one 



default observed 


in grade 


C , frequencies of obligors in grao 


'es given 


7 


50% 


75% 


90% 


95% 


99% 


99.9% 


PB 


0.52% 


0.73% 


0.95% 


1.10% 


1.43% 


1.85% 



From the given numbers of defaults in the different grades it becomes clear that a stand-alone 
treatment of grade B would yield still much higher values^ for the upper confidence bounds. The 
upper confidence bound 0.52% of the confidence region at level 50% is almost identical with the 
naive frequency based PD estimate 2/400 = 0.5% that could alternatively have been calculated 
for grade B in this example. 

For determining the confidence region at level 7 for the PD pc, by the same rationale as in 
Section [21 the grade C must be considered a stand-alone portfolio. According to the assumption 
made in the beginning of this section, one default occurred among the nc obligors in C. Hence 
we see that the confidence region for pc is the set of all admissible values of pc that satisfy the 
inequality 

1-7 < E(T)Pc(l-Pc)"^"' = {l-PcT'' +ncPc{l-Pcr^-'. (3.5) 

1=0 

For obligor frequencies as assumed in example 1)2. 3|) . Table IHl exhibits some maximum solutions^ 
Pc of (|3.5() for different confidence levels 7. 

*At level 99.9%, e.g., 2.78% would be the value of the upper confidence bound. 

^If we had assumed that two defaults occurred in grade B but no default was observed in grade C, then we 
would have obtained smaller upper bounds for pc than for p_g. As this is not a desirable effect, a possible - 
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Table 6: Upper confidence hound pc of pc as a function of the confidence 
level 7. No default observed in grade A, two defaults observed in grade B, one 



default observed 


in grade 


C , frequencies of obligors in grao 


'es given 


7 


50% 


75% 


90% 


95% 


99% 


99.9% 


PC 


0.56% 


0.90% 


1.29% 


1.57% 


2.19% 


3.04% 



So far, we have described how to generahze the methodology from Section [21 to the case where 
non-zero default frequencies have been recorded. In the following section we investigate the 
impact of non-zero default correlation on the PD estimates that are effected by applying the 
most prudent estimation methodology. 



4 Example: Correlated Default Events 



In this section, we describe the dependence of the default events with the one-factor probit 
modelf that was the starting point for developing the risk weight functions given in IBCBS 
(2004a|)^. First, we use the example from Section [2 and assume that no default at all was ob- 
served in the whole portfolio during the last period. In order to illustrate the effects of correlation, 
we apply the minimum value of the asset corr elation t hat ap pears in the Basel II corporate risk 
weight function. This minimum value is 12% ()BCBSl . l2004al . paragraph 272). Our model, how- 
ever, works with any other correlation assumption as well. Likewise, the most prudent estimation 
principle could potentially be applied to others than the Basel II type credit risk model as long 
as the inequalities can be solved for pA, Pb and pc, respectively. 

Under the assumptions of this section, the confidence region at le vel 7 for p 4 is repre sented as 
the set of all admissible values of pA that satisfy the inequality (cf. 
2.1.2 and 2.5.1 for the derivation) 



Bluhm et al 



2003, Sections 



1 - 7 < 



riA+riB+nc 



dy, 



(4.1) 



where ip and <I> stand for the standard normal density and standard normal distribution function, 
respectively. <&~^ denotes the inverse function of ^, and p is the asset correlation (here p is chosen 
as p = 12%). Similarly to (|2.2a|) . the right-hand side of inequality (|4.1j) tells us the one-period 
probability of not observing any default among ua + ?^_B + nc obligors with average PD pA- 



conservative - work-around could be to increment the number of defaults in grade C up to the point where pc 
would take on a greater value than pb- Nevertheless, in this case one would have to make sure that the applied 
rating system yields indeed a correct ranking of the oblig ors. 

® According to De Finetti's theorem (see, e.g.. fPurrettl Il99a . Theorem (6.8)), assuming one systematic factor 
only is not verv restrictive^ 

^See lCordvl t20031 andlBCBS ll2004H) for the background of the risk weight functions. In the case of non-zero 
realized default rates BalthaziMr i 20o3V uses the one-factor model for deriving confidence intervals of the PDs. 
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Solving® Equation (|4.1|) numerically^ for the frequencies as given in (|2.l'{j) leads to Table [3 with 
maximum solutions pA of (|4.1j) for different confidence levels 7. 

Table 7: Upper confidence hound pA of pA, Pb of pb and pc of pc as a 
function of the confidence level No defaults observed, frequencies of obligors 
in grades given by (|2.3() . Case of correlated default events. 



7 


50% 


75% 


90% 


95% 


99% 


99.9% 


PA 


0.15% 


0.40% 


0.86% 


1.31% 


2.65% 


5.29% 


PB 


0.17% 


0.45% 


0.96% 


1.45% 


2.92% 


5.77% 


PC 


0.37% 


0.92% 


1.89% 


2.78% 


5.30% 


9.84% 



Comparing the values from the first line of Table [7| with Table ^ shows that the impact of 
taking care of correlations is moderate for the low confidence levels 50% and 75%. The impact is 
much higher for the levels higher than 90% (for the confidence level 99.9% the bound is even six 
times larger) . This observation reflects the general fact that introducing unidirectional stochastic 
dependence in a sum of random variables entails a redistribution of probability mass from the 
center of the distribution towards its lower and upper limits. 

The formulae for the estimations of upper confidence bounds for ps and pc can be derived 
analogously to (|4.H) (in combination with (|2.4a|) and (|2.5aj) ). This yields the inequalities 

1 - . < jy,) (l - *(f:lfc^))"""° ,y (4.2a) 

and 

i-.<r^fe)(i-*(i^^^^i=^))"%., (4.2b) 



—00 



to be solved for ps and pc respectively. The numerical calculations with 1)4. 2a() and 1)4. 2b|) do 
not deliver additional qualitative insights. For the sake of completeness, however, the maximum 
solutions Pb of (|4.2a|) and pc of ()4.2b|l for different confidence levels 7 are listed in lines 2 and 
3 of Table [71 respectively. 

Second, we apply our correlated model to the example from Section |21 and assume that three 

defaults were observed during the last period. In analogy to Equations 1)3. 1|) . 1)3. 2|) and ()4.1() . the 

*See Appendix ^ Proposition I A. 21 for additional information. Taking into account correlations entails an 
increase in numerical complexity. Therefore, it might seem to be more efficient to deal with the correlation 
problem by choosing an appropriately enlarged confidence level in the independent default events approach as 
described in Sections |5| and |3 However, it remains open how a confidence level for the uncorrelated case, that 
"appropriately" adjusts for the correlations, can be derived. 

^The more intricate calcula tions for this paper were conducted by means of the software R (cf. 
R Development Core Teaiiii 2OO3I) . 
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confidence region at level 7 for pA is represented as the set of all values of pA that satisfy the 
inequality 

/oo 3 
^{y) 5] ("-^+7+"^) G{pA,p,yy (1 - (4.3a) 

where the function G is defined by 

G{p,p,y) = =^ . (4.3b) 

vi - P 

Solving (|4.3ap for pA with obligor frequencies as given in (|2.3p . and the respective modified 
equations for pB and pc yields the following results: 



Table 8: Upper confidence bound p a of pa, Pb ofpB andpc ofpc as a function 
of the confidence level No default observed in grade A, two defaults observed 
in grade B, one default observed in grade C , frequencies of obligors in grades 
given by H2.3() . Case of correlated default events. 



7 


50% 


75% 


90% 


95% 


99% 


99.9% 


PA 


0.72% 


1.42% 


2.50% 


3.42% 


5.88% 


10.08% 


PB 


0.81% 


1.59% 


2.77% 


3.77% 


6.43% 


10.92% 


PC 


0.84% 


1.76% 


3.19% 


4.41% 


7.68% 


13.14% 



Not surprisingly, the maximum solutions for PA, PB and pc increase if we introduce defaults in 
our example. Other than that, the results do not deliver essential additional insights. 



5 Potential Extension: Calibration by Scaling Factors 

One of the drawbacks of the most prudent estimation principle is that in the few defaults case, 
for all grades the upper confidence bound PD estimates are higher than the average default rate 
of the overall portfolio. This phenomenon is not surprising, given that we include all defaults of 
the overall portfolio in the upper confidence bound estimation even for the highest rating grade. 
However, these estimates might be regarded as too conservative by some practitioners. 

A potential remedy would be a scaling^'' of all of our estimates towards the central tendency 
(the average portfolio default rate) . We introduce a scaling factor K to our estimates such that 
the overall portfolio default rate is exactly met, i.e. 

^^"nt^C^n^"" ^ = ^^Portfolio. (5.1) 

The new, scaled PD estimates will then be 

PX,scaled = Kpx, X = A,B,C. (5.2) 

^"A similar scaling procedure has recently been suggested by Cathcart and BeniamirJ j2n 
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The results of the apphcation of such a scahng factor to our few defaults examples of Sections 
01 and are shown in Tables l9l and ITIH respectively. 

Table 9: Upper confidence bounds PA,scaied of pa, PB,scaied ofpB andpc,scaied of 
PC as a function of the confidence level 7 after scaling to the central tendency. 
No default observed in grade A, two defaults observed in grade B, one default 
observed in grade C , frequencies of obligors in grades given by (|2..S|) . Case of 
uncorrelated default events. 



7 


50% 


75% 


90% 


95% 


99% 


99.9% 


Central tendency 


0.375% 


0.375% 


0.375% 


0.375% 


0.375% 


0.375% 


K 


0.71 


0.48 


0.35 


0.30 


0.22 


0.17 


Pa, scaled 


0.33% 


0.31% 


0.29% 


0.29% 


0.28% 


0.27% 


P_B, scaled 


0.37% 


0.35% 


0.34% 


0.33% 


0.32% 


0.31% 


PCjScaled 


0.40% 


0.43% 


0.46% 


0.47% 


0.49% 


0.50% 



Table 10: Upper confidence bounds pA,scaied of pa, PB,scaied ofps andpc,scaied of 
PC as a function of the confidence level 7 after scaling to the central tendency. 
No default observed in grade A, two defaults observed in grade B, one default 
observed in grade C , frequencies of obligors in grades given by (|2.3|) . Case of 
correlated default events. 



7 


50% 


75% 


90% 


95% 


99% 


99.9% 


Central tendency 


0.375% 


0.375% 


0.375% 


0.375% 


0.375% 


0.375% 


K 


0.46 


0.23 


0.13 


0.09 


0.05 


0.03 


Pa, scaled 


0.33% 


0.33% 


0.32% 


0.32% 


0.32% 


0.32% 


P_B, scaled 


0.38% 


0.37% 


0.36% 


0.36% 


0.35% 


0.35% 


PC,scaled 


0.39% 


0.40% 


0.41% 


0.42% 


0.42% 


0.42% 



The average estimated portfolio PD will now fit exactly the overall portfolio central tendency. 
Thus, we loose all conservatism in our estimations. Given the poor default data base in typical 
applications of our methodology, this might be seen as a disadvantage rather than an advantage. 
By using the most prudent estimation principle to derive "relative" PDs before scaling them 
down to the final results, we preserve however the sole dependence of the PD estimates upon 
the borrower frequencies in the respective rating grades, as well as the monotony of the PDs. 

There remains the question of the appropriate confidence level for above calculation. Although 
the average estimated portfolio PD now always fits the overall portfolio default rate, the confi- 
dence level determines the "distribution" of that rate over the rating grades. In above example, 
though, the differences in distribution appear small, especially in the correlated case, such that 
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we would not explore this issue further. The confidence level could, in practice, be used to con- 
trol for the spread of PD estimates over the rating grades - the higher the confidence level, the 
higher the spread. 

However, above scaling only works if there is a non-zero number of defaults in the overall 
portfolio. Zero default portfolios would indeed be treated worse if we continue to apply our 
original proposal to them, compared to using scaled PDs for low default portfolios. 

A variant of above scaling proposal, that takes care of both issues, is the use of an upper 
confidence bound for the overall portfolio PD in lieu of the actual default rate. This upper 
confidence bound for the overall portfolio PD, incidently, equals the most prudent estimate for 
the highest rating grade. Then, the same scaling methodology as described above can be applied. 
The results of its application to the few defaults examples as in Tables El and ^1 are presented 
in Tables HU and ini 

Table 11: Upper confidence bounds pA,scaied of pa, PB,scaied ofps andpc,scaied of 
PC as a function of the confidence level 7 after scaling to the upper confidence 
bound of the overall portfolio PD. No default observed in grade A, two defaults 
observed in grade B, one default observed in grade C , frequencies of obligors 
in grades given by (|2.3|) . Case of uncorrelated default events. 



7 


50% 


75% 


90% 


95% 


99% 


99.9% 


Upper bound for portfolio PD 


0.46% 


0.65% 


0.83% 


0.97% 


1.25% 


1.62% 


K 


0.87 


0.83 


0.78 


0.77 


0.74 


0.71 


Pa, scaled 


0.40% 


0.54% 


0.65% 


0.74% 


0.92% 


1.16% 


PS, scaled 


0.45% 


0.61% 


0.74% 


0.84% 


1.06% 


1.32% 


PCjScaled 


0.49% 


0.75% 


1.01% 


1.22% 


1.62% 


2.17% 



Table 12: Upper confidence bounds pA,scaied ofpA, PB,scaied ofps andpc,scaied of 
PC as a function of the confidence level 7 after scaling to the upper confidence 
bound of the overall portfolio PD. No default observed in grade A, two defaults 
observed in grade B, one default observed in grade C , frequencies of obligors 
in grades given by (|2.3() . Case of correlated default events. 



7 


50% 


75% 


90% 


95% 


99% 


99.9% 


Upper bound for portfolio PD 


0.71% 


1.42% 


2.50% 


3.42% 


5.88% 


10.08% 


K 


0.89 


0.87 


0.86 


0.86 


0.86 


0.87 


Pa, scaled 


0.64% 


1.24% 


2.16% 


2.95% 


5.06% 


8.72% 


PB, scaled 


0.72% 


1.38% 


2.39% 


3.25% 


5.54% 


9.54% 


PC,scaled 


0.75% 


1.53% 


2.76% 


3.80% 


6.61% 


11.37% 
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As, in contrast to the situation of Tables IHl and [TOl in Tables ITTI and IT^ the overall default rate 
in the portfolio depends on the confidence level, we observe scaled PD estimates for the grades 
that increase with growing levels. Nevertheless, the scaled PD estimates for the better grades 
are still considerably lower than the corresponding unsealed estimates from Sections |21 and |^ 
respectively. For the sake of comparison, we provide in AnnexElthe analogous numerical results 
for the no default case. 

The advantage of this latter variant of the scaling approach is that the degree of conservatism 
is actively manageable by the appropriate choice of the confidence level for the estimation of 
the upper confidence bound of the overall portfolio PD. Moreover, it works for both the zero 
default and the few defaults case, and thus does not produce a structural break between both 
scenarios. Lastly, the results are less conservative than the ones of our original methodology. 

Consequently, we would propose to use the most prudent estimation principle to derive "relative" 
PDs over the rating grades, and subsequently scale them down according to the upper bound 
of the overall portfolio PD, which is once more determined by the most prudent estimation 
principle with an appropriate confidence level. 

6 Potential Extension: The multi-period case 

So far, we have only considered the situation where estimations are carried out on a one year 
(or one observation period) data sample. In case of a time series with data from several years, 
the PDs (per rating grade) for the single years could be estimated and could then be used 
for calculating weighted averages of the PDs in order to make more efficient use of the data. 
Proceeding this way, however, the interpretation of the estimates as upper confidence bounds at 
some pre-defined level would be lost. 

Alternatively, the data of all years could be pooled and tackled as in the one-year case. When 
assuming cross-sectional and intertemporal independence of the default events, the methodology 
as presented in Sections |21 and 13 can be applied to the data pool by replacing the one-year 
frequency of a grade with the sum of the frequencies of this grade over the years (analogous 
for the numbers of defaulted obligors). This way, the interpretation of the results as upper 
confidence bounds as well as the frequency-dependent degree of conservatism of the estimates 
will be preserved. 

However, when turning to the case of default events which are cross-sectionally and intertem- 
porally correlated, pooling does not allow for an adequate modelling. An example would be a 
portfolio of long-term loans, where in the intertemporal pool every obligor would appear several 
times. As a consequence, the dependence structure of the pool would have to be specified very 
carefully, as the structure of correlation over time and of cross-sectional correlation are likely to 
differ. 
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In this section, we present a multi-period extension of the cross-sectional one-factor correlation 
model that has been introduced in Section 01 We will take the perspective of an observer of a 
cohort of obligors over a fixed interval of time. The advantage of such a view arises from the 
possible conceptional separation of time and cross-section effects. Again, we do not present the 
methodology in full generality but rather introduce it by way of an example. 

As in Section^ we assume that, at the beginning of the observation period, we have got ua 
obligors in grade A, ub obligors in grade B, and nc obligors in grade C. In contrast to Section 
HJ the length of the observation period this time is T > 1. We consider only the obligors that 
were present at the beginning of the observation period. Any obligors entering the portfolio 
afterwards are neglected for the purpose of our estimation exercise. Nevertheless, the number of 
observed obligors may vary from year to year as soon as any defaults occur. 

As in the previous sections, we first consider the estimation of the PD pA for grade A. PD in 
this section denotes a long-term average one- year probability of default. Working again with 
the most prudent estimation principle, we assume t hat the PDs p a, Pb, and pc are equal, 
i.e. PA = Pb = PC = P- We assume, in the spirit of iGordvl (|2003), that a default of obligor 
i = 1, . . . , N = UA + ^^_B + nc in year t = 1, . . . , T is triggered if the change in value of their 
assets results in a value lower than some default threshold c as described below (Equation (|6.3|) ). 
Specifically, if Vi^t denotes the change in value of obligor z's assets, Vi^t is given by 



Vi,t = ^St + yA^Ci,t, (6.1) 

where p stands for the asset correlation as introduced in Section |1J St is the realisation of the 
systematic factor in year t, and ^j^t denotes the idiosyncratic component of the change in value. 
The cross-sectional dependence of the default events stems from the presence of the systematic 
factor St in all the obligors' change in value variables. Obligor i's default occurs in year t if 

Vi^i > c, . . . , V,,t-i > c, V,,t < c. (6.2) 

The probability 

P[y^,t <c] = Pi,t = P (6.3) 

is the parameter we are interested to estimate: It describes the long-term average one-year 
probability of default among the obligors that have not defaulted before. The indices i and t at 
Pi^t can be dropped because by the assumptions we are going to specify below pi^t will neither 
depend on i nor on t. To some extent, therefore, p may be considered a through-the- cycle PD. 

For the sake of computational feasibility, and in order to keep as close as possible to the Basel II 
risk weight model, we specify the factor varia bles St,t = 1, . . . , T , and = l,...,N,t = 

1, . . . ,r, as standard normally distributed (cf. iBluhm et al.l . 120031 1 . Moreover, we assume that 
the random vector (5i, . . . , St) and the random variables ^i^t, i = I, ■ ■ ■ ,N,t = 1, . . . , T, are 
independent. As a consequence, from (|6.1|1 follows that the change in value variables Vi^t are 
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all standard normally distributed. Therefore, implies that the default threshold^^ c is 

determined by 

c = <^-\p), (6.4) 
with $ denoting the standard normal distribution function. 

While the single components St of the vector of systematic factors generate the cross-sectional 
correlation of the default events at time t, their intertemporal correlation is effected by the 
dependence structure of the factors Si, ... , St- We further assume that not only the components 
but also the vector as a whole is normally distributed. Since the components of the vector are 
standardized, its joint distribution is completely determined by the correlation matrix 



(6.5a) 



1 r2,3 ••• r2,T 

rT-1,1 ■ ■ ■ rT-i,T-2 1 rT-i,T 
V rT,i ■ ■ ■ rT,T-2 rT,T-i 1 / 

Whereas the cross-sectional correlation within one year is constant for any pair of obligors, 
empirical observation indicates that the effect of intertemporal correlation becomes weaker with 
increasing distance in time. We express this distance-dependent behavior^'^ of correlations by 
setting in (|6.5aj) 

r,,t = s,t = l,...,r,s/t, (6.5b) 

for some appropriate < i9 < 1 to be specified below. 

Let us assume that within the T years observation period kA defaults were observed among the 
obligors that were initially graded A, ks defaults among the initially graded B obligors and kc 
defaults among the initially graded C obligors. For the estimation of pA according to the most 
prudent estimation principle, therefore we have to take into account k = kA + ks + kc defaults 
among obligors over T years. For any given confidence level 7, we have to determine the 
maximum value p of all the parameters p such that the inequality 

1 — 7 < P[No more than k defaults observed] (6.6) 

is satisfied - note that the right-hand side of (|6.6|) depends on the one-period probability of 
default p. In order to derive a formulation that is accessible to numerical calculation, we have 
to rewrite the right-hand side of 1)6. 6|) . 

^^At first sight, the fact that in our modei the default threshold is constant over time seems to imply that the 
model does not reflect the possibility of rating migrations. However, by construction of the model, the conditional 
default threshold at time t given the value Vi,t-i will in general differ from c. As we make use of the joint 
distribution of the V i ,. the refore rating migrations are implicitly taken into account. 

^ iBlochwitz et al.l (120041) proposed the specification of the intertemporal dependence structure according to 
Il6.5bll for the purpose of default probability estimation. 
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As the first step we develop an expression for obligor i's conditional probability to default 
during the observation period, given a realization of the systematic factors Si, . . . , St- From 
H6.1|) . (jnHJ, ()6.4() and by using of the conditional independence of the Vi^i, . . . ,Vi^T given the 
systematic factors, we obtain 

P[01igor i defaults I 5i,...,5t] = P [^jnin^ < | 5i, . . . , 5t] 

= 1 - F[^i,i > G{p, p, Si), ii,T > G{p, p,St)\Si,..., St] 

T 

= l-l[{l-^Gip,p,St))), (6.7) 
t=i 

where the function G is defined as in (|4.3b|) . By construction, in our model all the probabilities 
P[01igor i defaults | Si, . . . , St] are equal, so that, for any of the is, we can define 



7r(S'i, . . . , St) = P[01igor i defaults | Si, . . . , St] 

T 

= l-^{l-^{G{p,p,St))). 

t=i 

Using this abbreviation, we can write the right-hand side of 1)6. 6() as 

k 

P[No more than k defaults observed] = ^ E[P [Exactly t obligors default | 5i, . . . , St\\ 



(6. 



^=0 

^(^)E[7r(Si,...,ST)Ml-^(5i,---,ST)) 



N-e-i 



(6.8b) 



The expectations in ()6.8b|) are expectations with respect to the random vector {Si, . . . , St) and 
have to be calculated as T-dimensional integrals involving the density of the T-variate standard 
normal distribution with correlation matrix given by (|6.5aj) and ()6.5b|) . When solving (|6.6j) for p, 
we calculated the values of these T-dimensional integrals by means of Monte-Carlo simulation, 
taking advantage of the fact that the term 



k 

,N-i 



^ ( ^) vr(Si, . . . , St)' (1 - niSi, St))'"'' (6.8c) 



£=0 



can efficiently be evaluated by making use of (|A.1|) . 

In order to present some numerical results for an illustration of how the model works, we have 
to fix a time horizon T and values for the cros s-sectional cor relation p and the intertemporal 



correlation parameter i?. We choose T = 5 as iBCBSl (|2004aD requires the credit institutions 



to base their PD estimates on a time series with min imum length f ive years. For p, we choose 
p = 0.12 as in Section|ll i.e. again a value suggested bv lBCBSl (|2004a|). Our feeling is that default 
events with a five years time distance can be regarded as being nearly independent. Statistically, 
this statement might be interpreted as something like "the correlation of Si and S^ is less than 
1%". Setting i) = 0.3, we obtain corr[Si, S^] = i)^ = 0.81%. Thus, the choice •& = 0.3 seems to 
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be reasonable. Note that our choices of the parameters are purely exemplary, as to some extent 
choosing the values of the parameters is rather a matter of taste or of decisions depending on 
the available data or the purpose of the estimations. 

Table IT^ shows the results of the calculations for the case where no defaults at all were observed 
during five years in the whole portfolio. The results for all the three grades are summarized in 
one table. For arriving at these results, was first evaluated with = ua + nB + nc, then 
with N = ub + nc, and finally with = nc- In all three cases we set k = Q'm. ()6.8b|) in order to 
express that no defaults were observed. Not surprisingly, the calculated confidence bounds are 
much lower than those presented as in Table [7| demonstrating this way the potentially dramatic 
effect of exploiting longer observation periods. 

Table 13: Upper confidence bounds pA of pA, Pb of pb and pc of pc as a 
function of the confidence level 7. No defaults during 5 years observed, fre- 
quencies of obligors in grades given by 1)2. 3|) . Case of cross-sectionally and 
intertemporally correlated default events. 



7 


50% 


75% 


90% 


95% 


99% 


99.9% 


PA 


0.03% 


0.06% 


0.11% 


0.16% 


0.30% 


0.55% 


PB 


0.03% 


0.07% 


0.13% 


0.18% 


0.33% 


0.62% 


PC 


0.07% 


0.14% 


0.26% 


0.37% 


0.67% 


1.23% 



For Table [TH we did essentially the same computations as for Table IT^ the difference being that 
we assumed that during five years A;^ = defaults were observed in grade ^4, /cg = 2 defaults 
were observed in grade B, and fcc = 1 defaults were observed in grade C (as in Sections 
and |1] during one year). As a consequence, we had to set A; = 3 in ()6.8b|) for calculating the 
upper confidence bounds for pA and pb, as well as /c = 1 for the upper confidence bounds oipc- 
Comparing here with the results presented in Table El we observe again a very strong effect of 
taking into account a longer time series. 

Table 14: Upper confidence bounds pA of pA, Pb of pb and pc of pc as a 
function of the confidence level 7. During 5 years no default observed in grade 
A, two defaults observed in grade B, and one default observed in grade C . 
Frequencies of obligors in grades given by ()2.3() . Case of cross-sectionally and 
intertemporally correlated default events. 



7 


50% 


75% 


90% 


95% 


99% 


99.9% 


PA 


0.12% 


0.21% 


0.33% 


0.43% 


0.70% 


1.17% 


PB 


0.14% 


0.24% 


0.38% 


0.49% 


0.77% 


1.29% 


PC 


0.15% 


0.27% 


0.46% 


0.61% 


1.01% 


1.70% 
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7 Potential Applications 



The most prudent estimation methodology described in the previous sections can be used for 
a range of apphcations, both in a bank internal context as well as in a Basel II context. In 
the latter case, it might be of specific importance for portfolios where neither internal nor 
external default data are sufficient to meet the Basel requirements. A prime example might be 
Specialized Lending. In these high-volume, low-number and low-default portfolios, internal data 
are often insufficient for PD estimations per rating grade, and might indeed even be insufficient 
for central tendency estimations for the entire portfolio (across all rating grades). Moreover, 
mapping to external ratings - although explicitly allowed in the Basel context and widely used 
in bank internal applications - might be impossible due to the low number of externally rated 
exposures. 

The (conservative) principle of the most prudent estimation could potentially serve as an alterna- 
tive to the Basel slotting approach, subject to supervisory approval. In this context, the proposed 
methodology might be interpreted as a specific form of the Basel requirement of conservative 
estimations in case of data scarcity. 

In a wider, bank internal context, the methodology might be used for all sorts of low default 
portfolios. In particular, it could serve as a comple ment to other estimation method s, whether 



Schuermann and Hanson (12004) or others. 



this be mapping to external ratings, the proposals by 
As such, we see our proposed methodology as one additional source for PD calibrations, that 
should neither invalidate nor prejudge a bank's internal choice of calibration methodologies. 

However, we tend to believe that our proposed methodology should only be applied to whole 
rating systems and portfolios. The - at first sight imaginable - calibration of PDs of individual, 
low default rating grades by the most prudent estimation principle within an otherwise data 
rich portfolio seems infeasible because of the unavoidable structural break between average PDs 
(data rich rating grades) and upper PD bounds (low default rating grades). Similarly, we believe 
that the application of the methodology for back-testing or similar validation tools would not 
add much additional information, as the (e.g. purely expert based) average PDs per rating grade 
would normally be well below our proposed quantitative upper bounds. 



8 Open Issues 

For potential applications, a number of issues would need to be addressed. In the following, we 
list the ones that seem to be the most important to us: 

• Which confidence levels are appropriate? The proposed most prudent estimate could serve 
as a conservative proxy for average PDs. In determining the confidence level, the impact of 
a potential underestimation of these average PDs should be taken into account. One might 
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think that the transformation of average PDs into some kind of "stress" PDs, as done in 
the Basel 11 and many other credit risk models, could justify rather low confidence levels 
for the PD estimation in the first place (i.e. using the models as providers of additional 
buffers against uncertainty). However, this conclusion would be misleading, as it mixes 
two different types of "stresses" : the Basel II model "stress" of the single systematic factor 
over time, and the estimation uncertainty "stress" of the PD estimations. 

Nevertheless, we would argue for moderate confidence levels when applying the most pru- 
dent estimation principle, but according to another reasoning: The most common alterna- 
tive to our methodology, namely deriving PDs from averages of historical default rates per 
rating grade, yields a comparable probability of underestimating the true PD. As such, 
high confidence levels in our methodology would be hard to justify. 

• At which number of defaults should one deviate from our methodology and use "normal" 
average PD estimation methods (at least for the overall portfolio central tendency)? Can 
this critical number be analytically determined? 

• If the relative number of defaults in one of the better ratings grades is significantly higher 
than those in lower rating grades (and within low default portfolios, this might happen 
with only one or two additional defaults), then our PD estimates can turn out to be non- 
monotone. In which cases should this be taken as an indication for the non-correctness of 
the ordinal ranking? Certainly, monotony or non-monotony of our upper PD bounds do 
not immediately imply that the average PDs are monotone or non-monotone. Under which 
conditions would there be statistical evidence of a violation of the monotony requirement 
for the PDs? 

Currently, we do not have definite solutions to above issues. We believe, though, that some 
of them will involve a certain amount of expert judgment rather than analytical solutions. In 
particular, that might be the case with the first item. If our proposed approach would be used in 
a supervisory - say Basel 11 - context, supervisors might want to think about suitable confidence 
levels that should be consistently applied. 

9 Conclusions 

In this article, we have introduced a methodology for estimating probabilities of default in low 
or no default portfolios. The methodology is based on upper confidence intervals by use of the 
most prudent estimation. Our methodology uses all available quantitative information. In the 
extreme case of no defaults in the entire portfolio, this information consists solely of the absolute 
numbers of counter-parties per rating grade. 

The lack of defaults in the entire portfolio prevents reliable quantitative statements on both 
the absolute level of average PDs per rating grade as well as on the relative risk increase from 
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rating grade to rating grade. Within the most prudent estimation methodology, we do not use 
such information. The only additional assumption used is the ordinal ranking of the borrowers, 
which is assumed to be correct. 

Our PD estimates might seem rather high at first sight. However, given the amount of infor- 
mation that is actually available, the results do not appear out of range. We believe that the 
choice of moderate confidence levels is appropriate within most applications. The results can be 
scaled to any appropriate central tendency. Additionally, the multi-year context as described in 
Section El might provide further insight. 

Acknowledgment. The authors thank Til Schuermann, Claudia Sand, and two anonymous 
referees for providing useful hints on an earlier draft of this paper. 
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A Appendix 

This appendix provides some additional information on the analytical and numerical solutions 
of Equations and (gUJ). 

Analytical solution of Equation (|3.2|) . If X is a binomially distributed random variable with 
size parameter n and success probability then for any integer < < n we have 

T.r^)pHl- pT-' = nx<k] = i- v[Y <p] = ] (A.i) 

i=o Jo * (1 -tT " ^ dt 

with Y d enoting a beta d istributed random variable with parameters a = k + 1 and (3 = n — k 



Hinderer 



198C, Lemma 11.2). The beta distribution function and its inverse function 



(see, e.g., 

are available in standard numerical tools, e.g. in Excel. 

Direct numerical solution of Equation (|3.2() . The following proposition shows the existence 
and uniqueness of the solution of (|3.2j) . and, at the same time, provides us with initial values 
for the numerical root-finding (see (|A.2c() '). 

Proposition A.I Let < k < n be integers, and define the function fn,k '■ (0, 1) — > M 6y 

k 

fn,k{p) = E ( " ) ^^^1 - P)""'' P e (0' (A-2a) 

i=0 

Fix some < v < 1. Then the equation 

fnAP) = V (A.2b) 
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has exactly one solution < p = p{v) < 1. Moreover, this solution p{v) satisfies the inequalities 



1- \/v < p{v) < a/I - V. (A.2c) 
Proof. A straight-forward calculation yields 

-{n-k)rk)p'a-Pr-'-'- (A.3) 



dfn,k{p) _ / l„\/n\„kn „\n-k-l 



dp 

Hence fn,k is strictly decreasing. This implies uniqueness of the solution of HA.2b|) . The inequal- 
ities 

fnfiiP) < fn,k{p) < fn,n-l{p) (A.4) 

imply existence of a solution of HA.2b|) and the inequalities (|A.2c|) . □ 

Numerical solution of Equation (|4.1|) . For (|4.1|) we can derive a result similar to Proposition 
lA.ll However, there does not exist an obvious upper bound to the solution p{v) of ()A.5b|l as in 

Proposition A. 2 For any probability < p < 1, any correlation < p < 1 and any real number 
y define 

V ^ — P 

where we make use of the same notations as for Equation (|4.1() . Fix a value < u < 1 and a 
positive integer n. Then the equation 



V 



oo 

^{y){l-Fp{p,y)Ydy, (A.5b) 



with if denoting the standard normal density, has exactly one solution < p = p{v) < 1. This 
solution p{v) satisfies the inequality 

p{v) > 1 - (A.5c) 



Table 15: Upper confidence bounds pA,scaied ofpA, PB,scaied ofps and 

PC, scaled 

of PC as a function of the confidence level ^ after scaling to the upper confidence 
bound of the overall portfolio PD. No default observed, frequencies of obligors 
in grades given by (|2.3|) . Case of uncorrelated default events. 



7 


50% 


75% 


90% 


95% 


99% 


99.9% 


Upper bound for portfolio PD 


0.09% 


0.17% 


0.29% 


0.37% 


0.57% 


0.86% 


K 


0.61 


0.66 


0.60 


0.58 


0.59 


0.59 


Pa, scaled 


0.05% 


0.11% 


0.17% 


0.22% 


0.33% 


0.51% 


PB, scaled 


0.06% 


0.13% 


0.20% 


0.25% 


0.39% 


0.58% 


PCjScaled 


0.14% 


0.24% 


0.45% 


0.58% 


0.89% 


1.35% 
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Proof of Proposition lA.lL Note that - for fixed p and y - the function Fp{p,y) is strictly 
increasing and continuous in p. Moreover, we have 



= Wm. Fp{p,y) and 1 = Xim Fp{p^y). 

Equation (|A.6|) implies existence and uniqueness of the solution of HA.5b|) . 
Define the random variable Z by 

Z = Fp{p,Y), 



(A.6) 



(A.7) 



where Y denotes a stan dard nor i nallv distributed random variable. Then Z has the well-known 
Vasicek distribution fcf. IVasiceld . Il997^ ■ and in particular we have 

E[Z] = p. (A.8) 

Using (|A.7|) . Equation ()A.5b|) can be rewritten as 

V = E[(l-Zf]. (A.9) 

Since y ^ {1 — y)" is convex for < y < 1, by HA.8|) Jensen's inequality implies 

V = E[(l-Zf] > {1-pT. (A.IO) 

As the right-hand side of ()A.5b|) is decreasing in p, (jA.Scj) now follows from (|A.10() . □ 



B Appendix 



This appendix provides additional numerical results for the "scaling" extension of the most 
prudent estimation principle according to Section [5] in the case of no default portfolios. In the 
examples presented in Tables El find 1161 the confidence level for deriving the upper confidence 
bound for the overall portfolio PD, and the confidence levels for the most prudent estimates of 



Table 16: Upper confidence bounds pA,scaied of pa, PB,scaied ofps and pc,scaied of 
PC as a function of the confidence level 7 after scaling to the upper confidence 
bound of the overall portfolio PD. No default observed, frequencies of obligors 
in grades given by ()2.3() . Case of correlated default events. 



7 


50% 


75% 


90% 


95% 


99% 


99.9% 


Upper bound for portfolio PD 


0.15% 


0.40% 


0.86% 


1.31% 


2.65% 


5.29% 


K 


0.62 


0.65 


0.66 


0.68 


0.70 


0.73 


Pa, scaled 


0.09% 


0.26% 


0.57% 


0.89% 


1.86% 


3.87% 


PB, scaled 


0.11% 


0.29% 


0.64% 


0.98% 


2.05% 


4.22% 


PCjScaled 


0.23% 


0.59% 


1.25% 


1.89% 


3.72% 


7.19% 
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PDs per rating grade have always been set equal. Moreover, our methodology always provides 
equality between the upper bound of the overall portfolio PD and the most prudent estimate 
for PA according to the respective examples of Sections |H and |^ 
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